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Abstract 

In the presented study we investigated the second order endpoints of the lines of first order phase transitions which 
emerge for the QCD in the heavy and hght quark mass regime and for the three-dimensional three state Potts 
model with an external field. We located the endpoints with Binder cumulants and constructed the energy-like 
and ordering field like observables. The joint probability distributions of these scaling fields and the values of the 
Binder cumulant confirm that all three endpoints belong to the universality class of the 3-dimensional Ising model. 



1. Introduction 

The spontaneous breaking of global symmetries 
at low and their restoration at high temperature 
are common features of quantum field theories and 
statistical models. It is well known that universal 
properties at finite temperature phase transitions 
in (3+1) dimensional gauge theories are related to 
those in 3-dimensional spin models [1] . 

In the presented study we investigated the pro- 
posed connection between the SU(3) gauge theory 
with two heavy quark flavors and the three state 
Potts model with an external field [2] and in addi- 
tion the physically more relevant case with three 
light quark flavors [3] . In all three cases a first order 
phase transition line emerges and ends in a second 
order endpoint [4] . 

To determine the universality class of these end- 
points its precise location is mandatory. Moreover, 
observables which are sensitive to the universality 
classes are needed. They are constructed using a 
method originally proposed to study the liquid gas 



transition point [5] . We extended this approach to 
the more complex case of a quantum field theory. 



2. The Potts Model 

The 3-dimensional three state Potts model with 
an external field is defined with the spin variables 
a, G {1,2, 3} by 

H = -/3^JK,a,)-/i^<5(a„ag) (1) 
EE -PE - hM 



(2) 



A non vanishing h > favors the magnetization in 
the direction of the spin ag . The partition function 
on a finite lattice of size is then given by 



Z(/3,/i,L) = ^exp(-H). 



(3) 



The pseudo-critical line [3pc,L{h) is defined by the 
maxima of the susceptibility 
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^^ = -{{M^)-{M)^) 

A finite size analysis of xl leads to he G 
[0.0005,0.0010] [2] which is not sufficient for a 
precise location of the endpoint. 



3. Quantum Chromo Dynamics 

To simulate the QCD we used the standard stag- 
gered fermion action 



{L) = N-'{J2Tt n t^x(n,n4) 



(10) 



n 714=1 



in the limit iriq — > oo for the deconfinement tran- 
sition and the chiral condensate 

n/ 



(Trif-i) 



(11) 



in the limit —>■ for the chiral transition. 



SF{^,i',U,mq) = J2'^(^)^n^{U,mg)^{m) (4) 4. The Simulation Parameters 



with the fermion fields ip, ip and the link vari- 
ables t/^(n) = e-^9'^^''(")e5t/(3) in the staggered 
fermion matrix 

K^niU^mq) = 5nmmq + ^(_l)'U + ...+nM-i 

X \ {^n+p.,mU^{n) - dn,m+ixUl{m)) . (5) 

For the gauge action we choose the standard Wil- 
son plaquette action 

Sg{U)=Y. ^ jl-lReTrC/^.j. (6) 

Integrating out the fermion fields we obtain the 
partition function on a N'^ x A^^- lattice 



Z= I \ \ dU^{n) e 

n,iJ, 

with 



_gQCL 



(7) 



SQCD ^ |_/J5g(f/) + !yTVlnA:^^(?7,mg)} . (8) 

and n/ as the number of quark flavors. Expanding 
Trln a: = -Y.it TrM' with K = I - kM and 
K = l/2TOg one can show that for heavy quarks 
and iV^=4 

y n 



(9) 



i.e. the quark mass enters the action with 
The order parameters of the QCD phase transi- 
tions are the Polyakov loop 



For the Potts model we used the Wolff cluster 
algorithm with a ghost spin in order to implement 
the external field and to simulate on lattices of size 
y = 40^ — 70"^. We define a new configuration after 
1000 configurations to get autocorrelations of the 
energy E between 7 and 25 close to the pscudo 
critical points. We simulate at 3 to 4 different /3- 
values for every h and perform 10000 updates at 
each (/3, /i)-pair of couplings. In the analysis of the 
data we make use of the Multihistogram-Method 
[7] to interpolate between /? and h values. 

The QCD simulations are performed with the 
Hybrid R algorithm. We used lattices of size 
8^ X 4, 12''' X 4, 16'^ X 4 and for the heavy quark 
regime an additional size of 24^ x 4. We simu- 
lated with three mass degenerate quark flavors of 
masses m,a = 0.03, 0.0325, 0.035, 0.04 close to the 
chiral and two mass degenerate quark flavors of 
niqa = 1.5, 1.7, 1.8, 2.0 close to the deconfinement 
transition. For 3-4 different /?- values per mass pa- 
rameter we perform (1 — 3)/(3 — 7) x 10^ updates 
at small/large bare quark masses. In the analysis 
a reweighting in the /3-direction has been used. 



5. Constructing the Scaling Fields 

In order to construct observables for the Potts 
model which are more sensitive to the universality 
class of the endpoint we make the assumption 



M = M + sE, £ = E + rM 



(12) 



2 



and rewrite the Hamiltonian in terms of these new 
fields 

H = -t£- (13) 
where the new couphngs are given by 

^=^—{h-r(3), T = -^{p-sh). {U) 



1 



1 



Identifying r with the direction of the pseudocrit- 
ical Hne at the endpoint and demanding that the 

cncrgy-likc fluctuations and those of the ordering 
field are uncorrelated we fix the parameter r and s 

by 



dpcih) 



dh 



(SMSS) = 



(15) 



h=hc 



where 6X = X — {X). In QCD we define analo- 
gously for the light quarks 



M = il>ip + sSg, 
and similarly 

M = L + sSg, 



£ = Sg + riptp 



£ = SG+rL 



(16) 



(17) 



for heavy quark masses. Rewriting the action 
gQCD jj-^ terms of these variables is not possible, 
but the new energy-like and ordering field-like 
variables £ and M. define an effective Hamiltonian 



n 



t£ + ^M 



(18) 



which will control the universal behavior in the 
vicinity of the critical points in QCD. We similarly 
make a linear ansatz for the new couplings 



T={p-p,)+A{h-hc), 

^ = {h-hc) + B{p-l3,) 



(19) 
(20) 



with h = m{m ^) for light (heavy) quark masses. 
We use the constraints 



= 

and 
B = 



d{M) 



We 



dT 



= {SM S£)sQCD 



(21) 



dp 



(22) 



to fix three of the four constants r,s and A,B. 



Using a scaling ansatz n, = + -p- we found 
for the Potts model r^o = -0.690(2) « -s. Thus 
the transformation is just a rotation. In the case of 
the QCD a size dependence of r and s is observed 
but a quantitative scaling analysis was not feasible. 
We therefore used results from the largest lattice. 
Close to the endpoint we define for the chiral region 
(s, r) = (0.41 ± 0.52, 0.550 ± 0.007) and for the de- 
confinemt transition (s, r) = (0.67±2.48, -0.223± 
0.043). 



6. Locating the Endpoints 

For a precise location of the critical points we 
use the Binder cumulants 



Bs,L 



B, 



4,L 



(23) 



The cumulant B^^l should vanish for (3=(3pc and 
the cumulant B^^ l should take on a volume inde- 
pendent value ,64 c at the endpoint. This value is 
1.602(2), 1.242(2) or 1.092(3) for the 3-dimensional 
Ising, 0(2) or 0(4) universality classes, respec- 
tively. We found for the Potts model [2] a value of 
1.609(4)(10) at (0.54938(2),0.000775(10)) where 
the first error is a statistical error and the sec- 
ond one is due to uncertainties in locating the 
endpoint. For the QCD we found [3] a value of 
1.547(46) at (5.6756(5), 1.74(4)) and 1.639(24) at 
(5.1458(5),0.033(1)) in the heavy and Hght quark 
mass regimes respectively (Fig.l). This strongly 
suggest that all three endpoints fall in the Ising 
universality class. 



7. Joint Probability Distributions 

With the finite size scaling ansatz [6] for joint 
probability distributions of £ and M. 

Pl{£,M)xP£m {AjS£,A+^6M,A£T,AMC) (24) 

with A£ = a£L^/'' ,AM=aML'^~'^^'' and AmAXi = 
A£Ag=V one obtains the universal function 



3 



Pl{£, M) = A+A+,KM (A+<5£, A+^6M) (25) 

at ^ = T = (Fig. 2). Comparing these distributions 
with those obtained for other universahty classes 
we again conclude that the endpoint belongs to the 
Ising universality class. 



References 



1.70 
1.68 
1.66 
1.64 
1.62 
1.60 
1.58 
1.56 
1.54 f 
1.52 

0.00072 0.00074 0.00076 0.00078 0.0008 0.00082 



70^^ ■ 






* 50, 
40^ 


(a) 


^^^^^^^^^^^^^^^^ 






Ising 


*"^,A'''^X-^^ 
















h 



2.4 




24^4 ' — ■ — ' 






B4 


16^x4 (b) 




2.2 




12^4 ' — * — ' T 










2.0 
















1.8 

1.6 






- t- ■ 








Ising 


1.4 








1.2 






-4 - 

m 




[1] 

[2] 
[3] 
[4] 

[5] 
[6] 

[7] 



B. Svetitsky, L. Yaflfe, Nucl. Phys. B210 (1982) 423. 
F. Karsch, S. Stickan, Phys. Lett. B488 (2000) 319. 



F. Karsch, E. Laermann, Ch. Schmid, ^iep-lat/0107020 

S. Gavin, A. Gocksch, R.D. Pisarski, Phys. Rev. D 49 
(1994) 49. 

J.J. Rehr, N.D. Mermin, Phys. Rev. A 8 (1973) 472. 

N.B. Wilding, A.D. Bruce, J. Phys. Cond. Mat. 4 
(1992) 3087. 

A. M. Ferrenberg, R.H. Swendsen, Phys. Rev. Lett. 61 
(1988) 2635. 



0.030 



0.035 



0.040 



0.045 




Fig. 1. i34,i of the Potts model (a) and for the QCD with 
heavy (b) resp. light (c) quark masses. (Lines are drawn to 
guide the eyes) 



Fig. 2. Histogram "pcM. of the Potts model (a) at the 
endpoint and of the QCD with m = 1.8 (b), m = 0.325 (c) 
close to the end points. 
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